NASA Technical Memorandum 100086 


Numerical Study of Unsteady 
Viscous Hypersonic Blunt Body 
Flows with an Impinging Shock 

G. H. Klopfer, H. C. Yee and P. Kutler 


(NASA-TH- 100096) NUMERICAL STUDY OF N88-22650 

UNSTEADY VISCOUS HYPERSONIC BLUNT BODY FLOWS 
WITH AN INPINGING SHOCK (NASA) 10 p 

CSCL 12A Unclas 

G3/64 0142322 


April 1988 


NASA 

National Aeronautics and 
Space Administration 


NASA Technical Memorandum 100096 


Numerical Study of Unsteady 
Viscous Hypersonic Blunt Body 
Flows with an Impinging Shock 

G. H. Klopfer, NEAR Inc., Mountain View, California 

H. C. Yee 

P. Kutler, Ames Research Center, Moffett Field, California 


April 1988 


ru/\sA 

National Aeronautics and 
Space Administration 

Ames Research Center 

Moffett Reid, California 94035 




Preprint for the Proceedings of the Eleventh International Conference on Numerical Methods in Fluid 
Dynamics, June 27 - July 1, 1988, Williamsburg, Virginia. 


NUMERICAL STUDY OF UNSTEADY VISCOUS HYPERSONIC 
BLUNT BODY FLOWS WITH AN IMPINGING SHOCK 


G.H. Klopfer 1 

NEAR, Inc., Mountain View, CA 94043 
H.C. Yee 2 and P. Kutler 3 

NASA Ames Research Center, Moffett Field, CA 94035 


Abstract 

A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated 
by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of 
a non- stationary oblique shock impinging on the bow shock of a blunt body. Studies indicate that the 
unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases 
investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement 
locations than for the steady flow cases. 


Introduction 

The recent development of high resolution total variation diminishing (TVD) schemes [1-3] has made it 
possible to accurately simulate the unsteady Navier- Stokes equations for complicated shock interactions at 
hypersonic Mach numbers. Such flows are of current interest in hypersonic aerodynamics. In a previous 
study [4] the steady viscous hypersonic blunt body flows with impinging shock waves were investigated. 
Some numerical simulations by other numerical methods for steady blunt body flows with impinging shocks 
have been reported in [5-8]. Very little work has been done for viscous unsteady shock interaction flows at 
hypersonic speeds, in particular, numerical simulations by implicit methods. The objective of this paper 
is to study the unsteady viscous hypersonic blunt body flows with non- stationary impinging shocks. The 
rapidly changing shock waves and shear layers that occur and interact with each other make this a difficult 
problem to simulate by classical numerical methods. Here a two-dimensional code consisting of an implicit 
second-order accurate (time and space) TVD-type algorithm is used to simulate the flow field. The scheme 
is based on an existing TVD algorithm [9,10] for transonic and supersonic flows which has been extended 
recently to hypersonic and equilibrium real gas flows [2,3]. The specific objective of the study is to determine 
the unsteady shock interference patterns and to compare them with the corresponding patterns for steady 
flows at various hypersonic Mach numbers. Another objective is to determine the transient surface pressures. 


The Physical Problem 

The physical problem examined is the unsteady hypersonic viscous two-dimensional blunt body flow with 
an impinging shock. A schematic of the computational domain of the flow field is shown in figure L The 
blunt body has a thickness of D and a nose radius of The circumferential angle, 0, is measured from the 
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symmetry line as depicted in the figure. This physical model can represent many practical applications. One 
example is the bow shock of a hypersonic aerodynamic vehicle interacting with the shocks of the leading edge 
of a wing or the cowl lip of the engine inlet. The interaction causes complicated shock-on-shock patterns in 
steady flows. Such steady interactions were studied in [4]. However, even under steady cruise conditions the 
vehicle bow shock fluctuates about some mean position, and thus the interaction is unsteady. The thickness 
of the cowl lip or wing is thin compared to the shock layer thickness, and even a small shock fluctuation 
causes a large excursion of the impinging shock in the vicinity of the inlet cowl or wing. This unsteady 
interaction is modeled by a non-stationary impinging shock (vehicle bow shock) moving downward with 
constant speed across the bow shock of the blunt body which can represent the cowl lip or wing leading 
edge. For steady flows, various kinds of shock interactions occur depending on the freestream conditions and 
impingement shock angle and location. Based on experimental data at supersonic speeds, Edney [11] has 
identified six different kinds of interactions labeled Type I through Type VI. For the Type IV interaction 
for steady flows, the peak surface pressure is many times larger than for the blunt body flow case. One 
motivation for the study is to see if the same type of amplification of the peak surface pressure values also 
occurs for the unsteady interaction. 


Numerical Results 

Due to space limitations, numerical results obtained for the unsteady laminar thin layer Navier- Stokes 
flows of a perfect gas are only summarized here. At present there are no experimental data available for the 
unsteady interactions to validate the numerical simulations. Therefore the unsteady results are compared 
to each other at different Mach numbers and to the steady state computations of reference 4 at the same 
Mach numbers. 

Figure 2 presents the Mach contours at six instances of the diffraction process of an unsteady viscous 
computation. An impinging shock at an angle of 22.75° relative to the freestream moves down across the 
bow shock of the blunt body. The impingement shock velocity is 10% of the freestream velocity (Moo = 15). 
The freestream Reynolds number based on the diameter is Rep = 186,000 and the freestream temperature 
is Too = 255.6° -K - . The Mach contours ranging from 0 to 15 in increments of 0.1 are shown in figure 2. 
The indicated times of each frame are normalized with the freestream velocity and cowl lip thickness. The 
surface pressure distributions for the last three instances compared with the steady-state case are given in 
figure 3. The first three surface pressure distributions are not given here and are similar to the steady case 
[4]. The most critical condition in terms of surface pressures no longer seems to occur for the Type IV 
shock interaction, as it does for the steady state case. Experimental data [11] indicates that the rates can 
be as much as 1000% (10 times) higher than for the non- interfering blunt body flow. Note that the shock 
impingement locations of the unsteady and steady computations do not coincide. 

For the conditions used in this computation only the first five types of shock patterns identified by Edney 
are apparent in the sequence. Due to the unsteadiness, high Mach numbers, and high impingement shock 
angles the flow patterns are much more complex than those given by Edney. The corresponding Mach 
contours for the steady case at Moo = 15 are shown in figure 11 of reference 4. For both the steady and 
unsteady cases the same type of shock interference patterns occur although the unsteady sequence exhibits 
a much more complicated pattern within the shock layer. The separated boundary layer apparent for the 
unsteady Type I interaction does not. appear for the steady Type I case with the same shock impingement 
angle. Also the supersonic jets flowing parallel to the blunt body surface and below the impingement point 
in the Type III and IV unsteady flows do not occur for the steady state cases. However for the Type V shock 
interference, where a new 2-D blunt body flow forms with the freestream conditions of those of a supersonic 
wedge flow, there are very little differences between the steady and unsteady flows. 

The last set of results is shown in figures 4 and 5 for an unsteady Navier- Stokes computation at Mach 5.94. 
The freestream Reynolds number based on the diameter is Re d = 186,000 and the freestream temperature 
is Too = 59.5°ff. To demonstrate the differences at the two Mach numbers a similar sequence is used as for 
the Mach 15 case. For the Type I and II flows the shock patterns differ from the Mach 15 case for identical 
impingement shock location. However, for Type III to Type V flows, the shock patterns and locations are 
the same. Another notable difference at the two Mach numbers is that for the Type I flow, the transmitted 
shock causes the wall boundary layer to separate and thus becomes a lambda shock at the higher Mach 
number. The boundary layer does not separate at the lower Mach number. The boundary layer separation 
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has not been observed for the steady flows computed so far. The surface pressure distributions compared 
with the steady case are given in figure 5 for the last three of the six instances shown in figure 4. The 
Mach contours for the steady case can be found in figure 5 of reference 4. Similar to the Mach 15 case, the 
unsteady peak surface pressures no longer seem to occur for the Type IV shock interaction as they do for 
the steady flows. 

Comparison of the computed results for the surface pressures of the steady shock interaction with the 
experimental data of Keyes given in [12] can be found in reference 4 and is shown in figure 5. A complete 
description of the steady state results for various Mach numbers is also given in reference 4, where in 
particular the numerical computations are compared to the various experimental data in terms of the surface 
pressures and heat transfer rates. The comparison serves to validate the code used for the present study. 


Conclusions 

Numerical simulation of the unsteady Navier-Stokes equations with a high resolution time-accurate implicit 
TVD scheme for predicting the complicated shock-on-shock interaction on a blunt body with a non-stationary 
impinging shock has been carried out. It has been demonstrated that similar types of shock interactions 
occur for the unsteady flows as for the steady cases identified by Edney. However, for the sample cases 
studied in this paper the peak surface pressures do not seem to occur during the Type IV interaction as for 
the steady case. The significance of this fact is that the Type III interaction occurs for a much broader range 
of impingement shock angles and positions than for the Type IV interaction. Also the details within the 
shock layer differ appreciably. For example boundary layer separation and transient supersonic jets flowing 
parallel to the blunt body surface occur for the unsteady but not for the steady shock interactions. 
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Schematic of the computational 
domain. 



Fig. 2 Mach contours at six instances of the diffraction process with M <*, = 15 and Rep = 186,000. The 
indicated times (f*) of each frame are normalized with the freestream velocity and blunt body 
thickness. 
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Fig. 5 The surface pressure distributions for Types III, IV, and V shock interactions at M = 5.94. The 
unsteady sequence is for the last three instances shown in figure 4. The steady sequence is for the 
last three views shown in figure 5 of reference [4j. 
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